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Abstract: We study the impact of the cosmological parameters uncertainties on the measure- 
ments of primordial non-Gaussianity through the large-scale non-Gaussian halo bias effect. While 
this is not expected to be an issue for the standard ACDM model, it may not be the case for 
more general models that modify the large-scale shape of the power spectrum. We consider the 
so-called local non-Gaussianity model, parametrized by the /nl non-Gaussianity parameter which 
is zero for a Gaussian case, and make forecasts on /nl from planned surveys, alone and combined 
with a Planck CMB prior. In particular, we consider EUCLID- and LSST-like surveys and forecast 
the correlations among /nl and the running of the spectral index a^, the dark energy equation 
of state w, the effective sound speed of dark energy perturbations c^, the total mass of massive 
neutrinos Mi, — ^TOi^, and the number of extra relativistic degrees of freedom N^,^^ . Neglecting 
CMB information on /nl and scales k > 0.03/i/Mpc, we find that, if N'^J''^ is assumed to be known, 
the uncertainty on cosmological parameters increases the error on /nl by 10 to 30% depending 
on the survey. Thus the /nl constraint is remarkable robust to cosmological model uncertainties. 
On the other hand, if A^^°' is simultaneously constrained from the data, the /nl error increases by 
~ 80%. Finally, future surveys which provide a large sample of galaxies or galaxy clusters over 
a volume comparable to the Hubble volume can measure primordial non-Gaussianity of the local 
form with a marginalized l-tr error of the order A/nl ^2 — 5, after combination with CMB priors 
for the remaining cosmological parameters. These results are competitive with CMB bispectrum 
constraints achievable with an ideal CMB experiment. 
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1. Introduction 

Tests of deviations from Gaussian initial conditions offer an important window into the 
very early Universe and a powerful test for the mechanism which generated primordial 
perturbations. While standard single- field slow-roll models of inflation lead to small de- 
partures from Gaussianity, non-standard scenarios allow for a larger non-Gaussianity (NG) 
level (e.g. [|l|, ^, ^, and refs. therein). 

In particular, large NG can be produced if any of the conditions below is violated: a) 
single field, h) canonical kinetic energy c) slow roll and d) adiabatic (Bunch-Davies) initial 
vacuum state. The type of NG arising in standard inflation reads ^, ^, ^ 

^ = </' + /nl(0'-(0')) , (1.1) 

where $ denotes Bardeen's gauge-invariant potential, which, on sub-Hubble scales reduces 
to the usual Newtonian peculiar gravitational potential up to a minus sign, and (j) denotes a 
Gaussian random field. The NG parameter /nl is often considered to be a constant, yielding 
NG of the local type with a bispectrum which is maximized for squeezed configurations 
NG of the local type is generated in standard infiationary scenarios (where /nl is expected 
to be of the same order of the slow-roll parameters) as well as in multi-field inflationary 
scenarios 

The standard observables to constrain NG are the Cosmic Microwave Background 
(CMB) and the Large-Scale Structure (LSS) of the Universe. Traditionally, the most 



^Note that Eq. (1.1) is not general, i.e. there is a plethora of possible deviations from Gaussianity arising 



in the different inflationary scenarios proposed in the literature. 
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popular method to detect primordial NG has been to measure the bispectrum or the three- 
point correlation function of the CMB ^, 10|, while the LSS bispectrum has been shown 
to be sensitive to primordial NG only at high redshift 11, |l^, |l^, p!^ ]. 

Other powerful techniques to measure NG are based on weak lensing tomography |15| , 
Integrated Sachs- Wolfe effect (ISW) [0, abundance |l|, 0, 0, |2| and clustering 
1 23, 24 1 of rare events such as density peaks, since they trace the tail of the underlying 
distribution. 



Refs. |25| and |2^ (hereafter MV08) showed that primordial NG affects the clustering 
of dark matter halos inducing a scale-dependent large-scale bias. This effect, which goes 
under the name of non-Gaussian halo bias, is particularly promising, yielding already 
stringent constraints from existing data pT] , p^ ] . Forthcoming constraints on NG exploiting 
the non-Gaussian halo bias are expected to be similar to those achievable from an ideal 



CMB survey |Tq]. These predictions have been confirmed by N-body simulations |25, 28 



2|, |0|. 

Forecasts for /nl constraints from the halo-bias have been carried out so far assuming 
perfect knowledge of all other cosmological parameters. While for a ACDM model this is 
expected to be a reasonable assumption, for more general models one may expect /nl to 
be degenerate with other parameters and thus to have a larger marginal error. Here we 
study the degeneracies among the large-scale non-Gaussian halo bias (for NG of the local 
type) and the cosmological parameters which affect the large-scale halo power spectrum, 
focusing on dark energy perturbations, massive neutrinos, number of relativistic species, 
and running spectral index, which can produce large deviations of the underlying cosmology 
from the minimal ACDM scenario. The paper is organized as follows. In we briefly 
review the analytic expressions of the non-Gaussian halo power spectrum generalized to 
redshift dependent transfer functions. The redshift dependence is due to the presence of 
both dark energy perturbations and massive neutrinos. In §^ we summarize the Fisher 
matrix formalism applied to the observed halo power spectrum. In §Q we describe the 
assumed fiducial cosmology. Finally, in §^ and §|6| we discuss the results and draw our 
conclusions. 



2. Non- Gaussian halo bias 



Here we summarize the derivation of |26|, extending it to the case of a redshift dependent 
transfer function. In Fourier space, the filtered linear over-density 6r is related to the 
primordial potential $(k) by the Poisson equation: 

Sr{K z) = | ^(^m^'^)^'^V H(fc)$(k) = Mnik, zMk) , (2.1) 

where T{k, z) denotes the matter transfer function (which is redshift dependent in the 
presence of massive neutrinos and/or dark energy perturbations), Wii{k) is the Fourier 
transform of the top-hat function of width R, Hq and llm,o ^-re the current values of 
the Hubble constant and the total matter energy density respectively, and D{z) = (1 -|- 



- 2 - 



z)~'^g{z)/g{0) is the linear growth-factor of density fluctuations, normalized to D{0) = 1 
with g(z) being the growth suppression factor for non Einstein-de Sitter universes. 
In this context, the non-Gaussian halo power-spectrum takes the form 

Ph(A:, z) = b\^{z, M)Pss{k, z) [1 + 4/nl6l,/.(^, M)P{k, z)] . (2.2) 

Here bL^h{z, M) = q5c{z) / {a\jD{z)) is the Gaussian lagrangian bias^ for dark matter halos 
of mass M |Q, a\i is the mass variance linearly extrapolated to z = 0, 5c{z) = ^c{z) / D{z), 
being Ac(z) the linear overdensity for spherical collapse, which can be considered as a 
constant ^c{z) = 1.686, even in the presence of dark energy g is a factor extracted 
from N-body simulations ([^] and references therein) and ^^^(A;,^) is the power spectrum 
of 5r (as it will be shown later in Eq. (|4.8| )). Finally, /3(/c, z) is defined as 



p{k,z) = — - / dkiklMR{ki,z)P^{ki) di^iMR(Va,z 



-1 



, (2.3) 



where a = kf + k'^ -t- 2kikfi. Here we adopt the so-called CMB /nl normalization where 
Eq. ( |1 . 1| ) is intended to be deep in the matter-dominated era^. Consequently, the non- 
Gaussian halo Lagrangian bias reads'^ 

b^%{z,M)c^bL,h{^,M)[l + 2h^bL,h{^,M)f3{k,z)]. (2.4) 

Making the standard assumption that halos move coherently with dark matter, the Eulerian 
bias is 6^; = 1 + b^. The halo power spectrum given by Eq. (^^), is connected directly 
to the underlying dark matter power spectrum and can be recostructed from the galaxy 
power spectrum using different techniques (e.g. [^]). It provides important information 
on the growth of structure, which helps in constraining dark energy and neutrino masses 
together with primordial non-Gaussianities. We shall exploit information from both the 
shape and the amplitude of the NG halo power spectrum up to scales k < 0.03/i/Mpc, 
where /3{k, z) has a negligible dependence on the halo mass (recall that for local non- 
Gaussianity I3{k,z) oc l/k'^). In addition, Eqs. (p.2|)-( |2^ ) are valid only on scales much 
larger than the Lagrangian radius of the halo p6| , and have been tested and calibrated 
on N-body simulations only on scales k < 0.03/i/Mpc ||2^. We believe that the approach 
followed here is conservative, since in our parameter forecasts we do not consider scales 
k > 0.03/i/Mpc in the halo power spectrum, which could provide much tighter constraints 
on cosmology, via e.g. Baryonic Acoustic Oscillations (BAO). In what follows, we shall use 
Eq. ( p.2| ) to propagate errors of the non-Gaussian halo power spectrum into errors of the 
cosmological parameters via a Fisher matrix approach, as described below. 



^In general, the Gaussian halo bias may have a non trivial dependence on both the halo formation 
redshift Zf and the observation redshift Zo |3^, However, for objects that did not undergo recent 
mergers, zj 3> Zo, or in the case of rapid mergers Zf ~ Zo for 5^ S> <J%j, i.e. large M and/or high Zf, the 



bias is well approximated by Eq. (2.2) 

3 



See |g8| for the large-scale structure-normalized conversions. 
*A more accurate expression is given by b^'j^{z, M) = bL,h{z, M)y'^l + 4/NL&L,h(2, M)l3{k, z)) 
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3. Methodology 



In this paper we adopt the Fisher matrix formaUsm to make predictions of the cosmological 
parameter errors including the NG parameter /nl- The Fisher matrix is defined as the 
second derivative of the hkehhood surface about the maximum. As long as the posterior 
distribution for the parameters is well approximated by a multivariate Gaussian function, 
its elements are given by Q 



= ixr [C-^C,o.C-^C,^ 



(3.1) 



where C = S + N \s the total covariance which consists of signal S and noise terms. 
The commas in Eq. ( |3.1| ) denote derivatives with respect to the cosmological parameters 
within the assumed fiducial cosmology^. 

To derive realistic parameter forecasts, we consider future redshift surveys, as Euclid^- 
and LSST''-like galaxy surveys. While for BAO surveys an accurate redshift measurement 
is crucial [43, 44 1, for our purposes, a precise redshift extraction is not needed, as it will be 



explained after Eq. (3.3). For the EUCLID-like survey we will assume a redshift coverage 
of 0.5 < z < 2 and a sky area of /sky = 20000 deg^. For the LSST-like survey, we will 
consider 0.3 < z < 3.6 and /sky = 30000 deg^ for redshift and area coverages, respectively. 
In order to explore the cosmological parameters constraints from a given redshift survey, 
it is mandatory to specify the measurement uncertainties of the halo power spectrum. In 
general, the statistical error on the measurement of -Ph(^) at a given wavenumber bin is 



given by [^| 



APh 



n 2 



2(27r)2 



1 + 



1 



n 2 



(3.2) 



where nh is the mean number density of dark matter halos, V is the comoving survey 
volume of the galaxy survey, and n is the cosine of the angle between k and the line-of- 
sight direction. 

To our purposes it is adequate to perform an angular average over /i. Thus, our Fisher 
matrix for the large-scale structure data is given by 



pLSS 



'87r2 



dpo 



dpi3 



"-h-Ph(fc) 
flhPhik) + 1 



1 2 



(3.3) 



where pa represents the chosen set of cosmological parameters. 

We divide the surveys in redshift bins of width Az = 0.1 (larger than standard pho- 
tometric and spectroscopic redshift errors), and set /cmax to be 0.03/i/Mpc and /cmin to be 
greater than In/ AV^^^ , where AV is the volume of the redshift shell. Conservatively, we 



^In practice, it can happen that the choice of parameterization makes the posterior distribution shghtly 
non-Gaussian. However, even in this case, the error introduced by assuming Gaussianity in the posterior 
distribution can be considered as reasonably small, and therefore the Fisher matrix approach still holds as 
an excellent approximation for parameter forecasts. 

®http:/ /sci.esa.int/euclid 

''http: / / www.lsst .org/lsst /science/scientist_dark_energy 
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do not consider here that scales larger than /cmin can be used by cross-correlating different 
shells. 

The effect of NG alters the broad-band behavior of Phik) on very large scales, where 
Phik) is unaffected by the precision with which the radial positions of the galaxies are 
measured. Thus, we can treat photometric and spectroscopic surveys on the same footing. 
Moreover, the requirement of surveying a large volume of the universe and sampling highly 
biased galaxies to beat shot-noise, which is a key point for BAO surveys, is also a bonus 
for constraining primordial NG 

Note, moreover, that the NG correction of the halo bias is boosted by the Lagrangian 
Gaussian halo bias factor itself. 

In particular, for the value of fcmin used here, we find that i-h(^min) — Ph{k = 
0.2/i/Mpc), thus, the shot-noise requirement for BAO surveys of nP{k = 0.2h/Mpc) > 1, 
implies that for all scales of interest here, nP ^ 1. We have checked that our results do 
not change if we impose nP ~ 3 at all scales. 

We compute as well the CMB Fisher matrix to obtain forecasts for the Planck satel- 
lite^. We follow here the method of [^], considering the likelihood function for a realistic 
experiment with partial sky coverage, and noisy data 



■<BB 



hTTr^EE 1 r'TTr^EE _ onTEr^TE 

TT-tEE'^e 'r ^£ Z.y^i ^BB r, I fTT fEE j-BB 



-JfTTf 
y -I sky -I ; 




sky qTTqEE _ f^QTE'^2 ~^ •> sky qBB ■^yJskyJsky J sky f K'^-^J 

and compute its second derivatives to obtain the corresponding Fisher matrix 

In Eq. ( |3.4D Cf^ = Cf^ + being Cf^ the temperature and polarization power 

spectra (X, Y = {T, E, B}) and A/^ the noise bias. Finally, is the fraction of observed 
sky which can be different for the T-, and 5-modes. In Eq. ( p. 51 ) L = ln£, and pp 
denote the cosmological parameters of the assumed model and form the vector p whose 
fiducial value is given by p. 

Combining the Planck and redshift survey Fisher matrices {F^p = F^^^^+F^^^) we get 
the joint constraints for p. The 1~ct error on p^ marginalized over the other parameters 
is (j{pa) = a/ {F~^)aa, being F~^ the inverse of the Fisher matrix. We then consider 
joint constraints in a two-parameter subspace (marginalized over all other cosmological 
parameters) to study the covariance between /nl and the other cosmological parameters 
considered in this work. 

Furthermore, in order to quantify the level of degeneracy between different parameters 
and /nl, we estimate the so-called correlation coefficients, given by 



(3.6) 



^ www.rssd.esa.int /PLANCK 
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where pa denotes one of the model parameters. When the coefficient \r\ = 1, the two 
parameters are totally degenerate, while r = means they are uncorrelated. 



Table 1: Correlation coefBcients for M,y-cosmology 









h 






Us 


w 


c2 


as 




LSST 
LSST+Planck 


/nl 
/nl 


0.50 
-0.22 


0.31 
-0.19 


0.36 
0.36 


0.35 
-0.12 


-0.03 
0.14 


0.41 
-0.07 


-0.32 
0.10 


-0.31 
-0.13 


0.03 
0.34 


EUCLID 
EUCLID+Planck 


/nl 
/nl 


0.55 
-0.17 


0.25 
-0.05 


0.31 
0.24 


0.36 
-0.07 


-0.06 
0.06 


0.39 
-0.07 


-0.30 
0.09 


-0.24 
-0.05 


0.09 
0.16 



Table 2: Correlation coefficients for A^^° -cosmology 









h 






Us 


w 






^rcl 


LSST 
LSST+Planck 


/nl 
/nl 


-0.06 
0.24 


0.34 
-0.48 


0.12 
0.61 


-0.13 
-0.24 


0.12 
0.19 


0.34 
0.54 


-0.39 
0.04 


-0.44 
-0.08 


-0.64 
0.69 


EUCLID 
EUCLID+Planck 


/nl 
/nl 


-0.05 
0.34 


0.34 
-0.57 


0.08 
0.76 


-0.20 
-0.26 


0.05 
0.21 


0.11 
0.64 


-0.26 
0.08 


-0.45 
-0.05 


-0.69 
0.77 



4. Model parameters 

The Fisher matrix approach propagates errors of the observed P^, see Eq. (|3.2| ), into 
errors of the cosmological parameters which characterize the underlying fiducial cosmology. 



According to the latest observations (e.g. [51| and refs. therein), our fiducial ACDM 
cosmological parameters are: 0^,0^^ = 0.1358, O.^ oh'^ = 0.02267, h = 0.705, A'^{ko) = 
2.64 X 10~^, Hg = 0.96, Us = and /nl = 0. Here ^mfl and ^b,o B^ie the total matter and 
baryon present-day energy densities, respectively, in units of the critical energy density of 
the Universe, h is given by Hq = lOO/i km Mpc~^ , where Hq is the Hubble constant, 
A^(A:o) represents the dimensionless amplitude of the primordial curvature perturbations 
evaluated at a pivot scale /cq, Ug is the scalar spectral index of the primordial matter 
power spectrum, assumed to be a power-law, and Og is the running of the scalar spectral 
index. We will consider two different fiducial models matching this ACDM cosmology 
as follows. We adopt the same values for the 7 "base" parameters. We do not consider 
primordial gravitational waves and assume a flatness prior, fix = 0, as predicted by long- 
lasting inflation models, so that ^deo = 1 ~ ^m,0) where and i^deo are, respectively, 
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Table 3: /nl l-c^ errors for M„ cosmology 



fixed parameter LSST LSST+PLANCK EUCLID EUCLID+PLANCK 



non-marginalized 


1.65 


1.65 


2.79 




2.79 


marginalized 


4.52 


2.11 


8.86 




3.12 




3.91 


2.06 


7.41 




3.07 


h 


4.29 


2.07 


8.58 




3.11 




4.22 


1.97 


8.43 




3.02 




4.22 


2.10 


8.27 




3.11 


Us 


4.52 


2.09 


8.85 




3.11 


w 


4.12 


2.11 


8.16 




3.11 




4.27 


2.10 


8.44 




3.11 


as 


4.30 


2.10 


8.61 




3.11 


M, 


4.52 


1.99 


8.83 




3.08 




4.05 


2.09 


8.13 




3.10 


as, Ml, 


4.28 


1.94 


8.57 




3.07 




4.26 


1.97 


8.35 




3.06 


ci,M^,as 


4.02 


1.93 


8.03 




3.06 




Table 4: /nl 


1-cr errors for A''^°'-cosmolo 


gy 




fixed parameter 


LSST LSST+PLANCK 


EUCLID 


EUCLID+PLANCK 


non- marginalized 


1.46 


1.46 


2.56 




2.56 


marginalized 


5.08 


2.56 


10.15 




4.79 




5.07 


2.49 


10.13 




4.51 


h 


4.77 


2.25 


9.53 




3.93 




5.05 


2.03 


10.12 




3.13 




5.04 


2.49 


9.94 




4.62 


Us 


5.05 


2.51 


10.13 




4.68 


W 


4.78 


2.16 


10.09 




3.67 


(? 


4.69 


2.56 


9.79 




4.77 


as 


4.57 


2.55 


9.06 




4.78 


^rcl 


3.88 


1.86 


7.32 




3.06 


Us, Cg 


4.16 


2.55 


8.70 




4.77 


as,Nl^' 


3.76 


1.76 


7.15 




2.92 




3.57 


1.86 


7.15 




3.05 




3.43 


1.76 


6.95 




2.91 



the present-day energy densities associated to the spatial curvature and to the dark energy 
component of the Universe, in units of the critical density. 
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A Gaussian prior of 5% on the present-day Hubble's constant Hq = 
lOO/i km s^^ Mpc~^ is assumed, following the results of |52]. While this uncertainty is 



comparable to the one achieved by recent WMAP-7yr data^ in the determination of Hq 
for the ACDM model |51|, this information will improve, as expected, the parameter con- 
straints on models different from the minimal ACDM model, such as models with dark 
energy perturbations, massive neutrinos and non-vanishing running of the spectral index. 

We also consider dark energy to be a cosmic fluid with clustering properties on the 
Gpc scale, described by an equation of state that we assume to be constant 

Pde I / ,1 1 \ 

w = — ='w|fid, (4.1) 



Pd, 



where where pde and pde represent respectively the pressure and energy density of the 
dark energy fluid, and we assume a fiducial value = —0.9 which lies well within the 
current 95% C.L. limits on a constant dark energy equation of state parameter w. The dark 
energy fluid will be also described here by an effective sound speed Cg which parametrizes 
the transition between the smooth and clustered dark-energy regimes 



d^P^, (4,2) 

OPde 

whit fiducial value c^|fid = 0.9. Dark energy perturbations will arise only if the dark energy 
equation of state parameter w is different from —1. 

We assume the power spectrum of primordial curvature perturbations, P'ji{k), to be 

^Uk) - = AUko) (-] . (4.3) 



where /cq = 0.002/Mpc and A'^{ko)\M = 2.64 x lO^^ |8|. 

The matter energy density Qm,o includes the neutrino contribution when neutrinos are 
non-relativistic 

^m,0 = ^c,0 + ^^6,0 + ^u,0 , (4.4) 

where 0,,^^ is related to the sum of neutrino masses Mj^ = ^ m,y as 

and the neutrino mass eigenstates are assumed to have a degenerate spectrum, i.e. the 
three neutrinos have the same mass. 

We will constrain the following set of "baseline" parameters 

Pa = {^b,oh'^,h,0.cflh'^,A'^{ko),ns,w,Cs,as,fNL} ■ (4.6) 
From here we specify two cosmological models: 



^http:/ /lambda.gsfc. nasa.gov/ 
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A^^*^^-cosmology, where neutrinos are effectively massless but the the number of rel- 
ativistic species N^,^^ can deviate from the standard value N^'^^ = 3.04. In this case 
the fiducial value A^^J'^'Im = 3.04 is chosen, fixing = const = 0. A''^''^ is given by 
the energy density associated to radiation 

fi,,o = J^7,o (l + 0.2271iV;<=i) , (4.7) 

where il^^o = 2.469 x lO^^/i^^ is the present-day photon energy density parameter 
for Tcn,b = 2.725 K [||]. 

Mi,-cosmology, where A'^^^^ is fixed at the fiducial value and Mj, is allowed to vary. 



In this case, being still consistent with current data [50, 51 1 we choose a fiducial 
value M,y|fid = 0.3 eV. This choice is motivated by the fact that for taking two-sided 
numerical derivatives, the fiducial M^, must be non-zero. It is well known that the 
error from cosmological observations on the neutrino mass depends somewhat on the 
fiducial mass chosen; from we estimate that around a fiducial = the error 
on = would increase by less than 20%. As it will be clear from Sec. ^ the 
effect of this correction will be negligible on the /nl error estimate. In addition most 
of the signal to constrain Mi, from LSS surveys will come from smaller scales (not 
considered here). 

At the CMB level, if neutrinos are still relativistic at the decoupling epoch, z ~ 1090, 
i.e. if the mass of the heaviest neutrino specie is rrii, < 0.58 eV, massive neutrinos do not 
affect the CMB power spectra, except through the gravitational lensing effect |^], and, 
as a consequence, the dark energy equation of state w is not degenerate with the neutrino 
mass. However, the limit on the the sum of the neutrino masses degrades significantly 
when the dark energy equation of state is a function of redshift, if the amplitude of the 
galaxy spectrum is used for getting constraints on w and M^, since dark energy and massive 
neutrinos both affect the growth rate of structures. 

In this work, both massive neutrinos and clustering properties of the dark energy 
perturbations are considered. In this scenario, the growth function of the dark matter 
perturbations is scale-dependent, even at the linear level. The overall effect induces a 



redshift-dependent transfer function ||5^, and the power spectrum of the linear 



density field, smoothed on a sphere of radius R, takes the form 

Pssik, z) = %;^^^r'°V |(fc)r^(fc, z)D\z) ( A j , (4.8) 

where D{z) is the scale independent linear growth-factor defined in The redshift- 
dependent transfer function T{k, z) is directly extracted from CAMB^'^[57| at each redshift 



z, in order to compute the Fisher matrix, given by Eq. (p.3|), within each redshift bin. 

Our analysis exploit exclusively the linear matter power spectrum, since we restrict 
ourselves to scales k ^ 0.03/i/Mpc, where the details of the halo occupation distribution of 



"http:/ /camb.info/ 
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galaxies are irrelevant. We focus here on dark matter halos with mass ~ 

10i2_ioi3 (where 

the lower mass limit is relative to the highest redshift of the survey) . 

Notice as well that in our Fisher matrix analysis we do not add constraints on the /nl 
parameter from CMB experiments, so that our forecasts on /nl result exclusively from 
future redshift survey measurements of the dark matter halo power spectrum on scales 
k ^ 0.03/i/Mpc, i.e. without including information from BAOs which will further reduce 
the forecasted errors and residual degeneracies. For all the reasons explained above, the 
results presented in the next section should be considered conservative. 

5. Results 

In this Section we present the predicted 1-a marginalized errors of the /nl parameter, and 
the /nl covariance with the remaining cosmological parameters considered in our Fisher 
matrix analysis. We show forecasts both from LSST and EUCLID data only, as well as 
the expected errors after combining the results from these two experiments with Planck 
forecasted errors. 

First of all, when considering forecasts from the redshift surveys alone, we expect that 
the /nl parameter will be correlated with all cosmological parameters which affect the 
amplitude and shape of Ph{k) at scales k ^ 0.03/i/Mpc. Tables show the correlation 
r (see Eq. ( |3.6| )) among /nl and the other cosmological parameters pa- We expect w and 
/nl to be correlated. In fact, at scales k ^ 0.03/i/Mpc, an increase of w produces a trend 
on the matter power spectrum which is opposite to the one produced in Ph{k) by increasing 
/nl- So the effect of positively increasing /nl can be mimicked by a larger w (see Fig. |l]). 

Likewise, /nl is correlated with H.hflh'^, ^cflh'^, and M^,, since, on the scales consid- 
ered, the larger these parameters are, the smaller the matter power spectrum is, and this 
competes with the rise of Ph{k) due to an increasing value of /nl- 

On the other hand, /nl is negatively correlated with both the running of the scalar 
spectral index and the effective sound speed of dark energy perturbations c^. In fact, 
if either as or Cg increase, the matter power spectrum Pss is modified in a very similar 
way to the halo power spectrum Ph{k) for a larger /nl value. Consequently, the effect 
of a positively increasing /nl can be mimicked by decreasing either Og or c^. For what 
concerns the effective number of relativistic species, the ft^h~^l^^ correlation is more 
complicated, as it depends on several factors. Naively, one would expect a positively 
correlation between these two parameters, since an increase in the number of relativistic 
particles should suppress the matter power spectrum, which can be compensated by 
increasing /nl- However, since the redshift evolution of the halo bias is also modified in 
this situation, the correlation coefficient may change its sign. 

When the Planck Fisher matrix information is added to the survey Fisher matrix, 
all degeneracies are either resolved or drastically reduced. In some cases, the correlation 
coefficient r can even change sign, see Tabs. |l[]2|. This change in the behavior of r arises 
either due to the presence of dominant parameter degeneracies affecting the CMB spectrum, 
or because of marginalization of a high-dimension parameter space down to two variables. 
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Figure 1: The solid blue line represents the NG halo power spectrum for /nl = 100 and halo 
mass M = 10^^ Mq, while the dotted blue line is the corresponding Gaussian halo power spectrum 
(/nl ~ 0), at redshift z = 0. The dashed red line represents the matter power spectrum calculated 
for w — —0.9, while the red three-dot-dashed line is the matter power spectrum for w — —0.2. The 
matter power spectra are evaluated at z = and normalized to the same amplitude at A; = 0.002. 
Note that A{\n Ph{k)) / Aw and A(lnPft,(fc))/A/NL have opposite sign in the range of k of interest 
here: it is clear from Eq. (3.3) that the two parameters can compensate each other, i.e. they are 
correlated 



In particular, it is worth noting here that, while /nl and Nl'^^ are negatively correlated if 
only galaxy survey data are considered, they are positively correlated after adding Planck 
priors. In fact, since both N^^^^ and ^cflh'^ are strongly, positively correlated at the CMB 
level via the equality redshift Zeq, then, a positive correlation between ^Icfih'^ and /nl 
automatically turns into a correlation of the same sign between N^^^ and /nl- Despite the 
residual non-zero correlation coefficients, one should bear in mind that the marginalized 
/nl errors decrease by a factor > 2 when the CMB prior is added. 

From our analysis, we conclude that the effective number of relativistic species is the 
main parameter affecting the constraints on /nl- 

The 1-a errors of /nL; for the two fiducial cosmologies considered here, are shown in 
Tabs. ^ and |^. Let us notice that the marginalized errors are significantly larger than the 
non-marginalized ones when only the LSS surveys are used for the forecasts. Nonetheless, 
the marginalized errors become comparable in magnitude to the non-marginalized errors 
when Planck priors are added, since the CMB mitigates the intrinsic degeneracies between 
/nl and the other cosmological parameters at the LSS level. It is also worth noting that the 
non-marginalized 1-a errors of /nl presented in this work are larger than the corresponding 
errors presented in jlC]. In fact we now consider less highly biased halos, with a fiducial 



bias parameter more in-line with the expected one for (blue) EUCLID galaxies 59 1. It 



Figure 2: 2-parameter /nl-Pq joint contours for tiie fiducial model with extra relativistic degrees of 
freedom N^°^ as described in the text, obtained after combining LSST (upper panels) and EUCLID 
(lower panels) data with Planck priors. The blue dotted line, the red dashed line and the orange 
dot-dashed line represent the 68% C.L., 95.4% C.L. and 99.73% C.L., respectively. The black solid 
line shows the 1-parameter confidence level at 1-a. 



is important to note that the /nl effect on the halo bias is modulated by bL^h{z, M), which 
depends crucially on the selected halo and its merging history 61]. 

Tabs. show as well the effect of each cosmological parameter on the /nl forecasts. 
Obviously the parameters which have a larger impact on the /nl errors are the ones more 
degenerated with it, and can be directly inferred from Tabs. |l[]2|. Moreover, since with the 
inclusion of the parameters Cg, My and A^^^', we have considered cosmologies which 
deviate substantially from the minimal ACDM model, we fix pairs/triplets composed by 
these parameters to show how much deviations from a ACDM cosmology can affect the 
/nl constraints. For the N^^^ model cosmology there is an important impact on the /nl 
marginalized errors from Qcflh'^, h, w and, in particular, from N^'^^. In summary, if N^^^ 
is assumed to be fixed, the uncertainties on the other cosmological parameters increase 
the error on /nl only by 10 to 30%, depending on the survey. If N^^^ is considered as an 
extra parameter to be simultaneously constrained from the data then the uncertainty in 
the underlying cosmology increases the /nl error by ~ 80%. 

In Figs. |-| we show the 2-parameter projected 68% C.L., 95.4% C.L. and 99.73% C.L. 
contours in the /nl-Pq sub-space with pa = w^c^^Os, M,y, N^'^^, obtained after combining 
LSST and EUCLID data with Planck priors for the two fiducial models considered in this 
work. The black line shows the 1-parameter confidence level at 1-a. The orientation of 
the ellipses reflects the correlations among the parameters shown in Tabs. 



6. Conclusions 

Deviations from non-Gaussianity, usually parameterized by the parameter /nl, offer a 
powerful tool to identify the mechanism which generates the seeds for the structures we 
observe currently in our Universe. 

Here we study the impact of the uncertainties of the cosmological parameters on the 
/nl errors expected for the case of local non-Gaussianity for the large-scale non-Gaussian 
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Figure 3: The same as in Fig. 2, for the fiducial model with massive neutrinos of total mass 
Afi/|fid — 0,3 eV, as described in the text. 



halo bias effect. We forecast the correlations among /nl and the remaining cosmological 
parameters (including the running of the spectral index Og, and the dark energy parameters 
w and c^) within two possible cosmological models. The first model contains massive 
neutrinos (hypothesis robustly confirmed by neutrino oscillation data), where the total 
neutrino mass is a parameter to be constrained by the cosmological data. The second 
model assumes massless neutrinos (or neutrinos with a mass too small to be relevant for 
the cosmological observations considered here) and allows for extra relativistic degrees of 
freedom Nl^\ which could be induced by the presence of sterile neutrinos, non minimally 
coupled quintessence fields, or even by the violation of the spin statistics theorem in the 
neutrino sector. 

We follow here a conservative approach, assuming that /nl is constrained exclusively 
from the very large scale halo power spectrum (i.e. we neglect CMB information on /nl); 
and restrict ourselves to scales k ^ 0.03/i/Mpc, without exploiting information e.g., from 
BAOs, which will further reduce degeneracies and forecasted errors. We present first the 
Fisher matrix forecasts for /nl assuming EUCLID- and LSST-like surveys for the two 
model cosmologies considered here. Then, we add the Planck Fisher forecasts for the 
remaining cosmological parameters to study the impact on the /nl correlations. 

The combined errors on /nl do not change significantly in the presence of a dark 
energy equation of state, massive neutrinos, running of the spectral index, or clustering 
of dark energy perturbations, which are the parameters we have particularly focused on, 
since they are expected to affect the matter power spectrum on large scales, and represent 
the main deviations from a minimal ACDM model. However, the errors on /nl are highly 
affected in the presence of extra relativistic degrees of freedom iV^^i. We find that if iV^^' 
is assumed to be fixed, the effect of the uncertainties on the other cosmological parameters 
increases the error on /nl only by 10 to 30% depending on the survey. If A'^,^'^' is considered 
as a parameter to be simultaneously constrained from the data, then the uncertainty in the 
underlying cosmology increases the /nl error by ~ 80%. We thus conclude that, except 
for the effect of N^'^^, the halo-bias /nl constraints are remarkably robust to uncertainties 
in the underlying cosmology. 
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One important point to discuss is the effect of the (Gaussian) halo bias, as its value 
boosts the effect of /nl on the halo power spectrum shape, and in our analysis it has been 
assumed to be known. The (Gaussian) halo bias depends strongly on the type of halos 
selected by the survey -whether they correspond to extremely high and rare peaks in the 
initial fluctuation field-, and on their accretion history. Errors on /nl may be improved 
-at least in principle- by up to a factor of two by optimizing the choice of tracers. 

The bias factor itself will need to be estimated from the survey, at the same time as 
the other cosmological parameters; the signal comes from scales much smaller than those 
used here, where the NG effect on halo bias is completely negligible. We estimate that the 
error on the (Gaussian) halo bias will be of the same order (in %) as the error on the linear 



growth factor / as a function of redshift, which is forecasted to be <~ 10% [£2|, |6§]. Such 
a residual uncertainly will therefore increase the /nl errors reported here by at most 10%. 

Let us recall that the purpose of this work is to show the main correlations between 
/nl and the other cosmological parameters, and to understand if these degeneracies can 
degrade dramatically the /nl errors. We have shown that, after the combination with 
Planck constraints on parameters different from /nLj the degeneracies get mostly broken, 
independently on the particular cosmological parameter, even without adding information 
from smaller scales corresponding to A: > 0.03/i/Mpc. Therefore we conclude that the /nl 
constraints are very robust against underlying cosmology assumptions. 

Finally, future surveys which provide a large sample of galaxies or galaxy clusters over 
a volume comparable to the Hubble volume (LSST, EUCLID) will measure primordial 
non-Gaussianity of the local form with a marginalized 1-a error of the order A/nl ~ 2 — 5, 
after combination with CMB priors for the remaining cosmological parameters. These 
results are competitive with CMB bispectrum constraints achievable with an ideal CMB 



experiment A/nl ~few |^, 65]. 
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